A pseudo-unitary ensemble of random 
matrices, PT-symmetry and the Riemann 

Hypothesis 



Zafar Ahmed and Sudhir R. Jain 

Nuclear Physics Division, Van de Graaff Building, 
Bhabha Atomic Research Centre, Trombay, Mumbai 400 085, India 

Abstract 

An ensemble of 2 x 2 pseudo-Hermitian random matrices is con- 
structed that possesses real eigenvalues with level-spacing distribution 
exactly as for the Gaussian unitary ensemble found by Wigner. By 
a re-interpretation of Connes' spectral interpretation of the zeros of 
Riemann zeta function, we propose to enlarge the scope of search of 
the Hamiltonian connected with the celebrated Riemann Hypothesis 
by suggesting that the Hamiltonian could also be PT-symmetric (or 
pseudo-Hermitian) . 



PACS Nos : 05.45.+b, 03.65.Ge 



Riemann Hypothesis (RH) states that all the nontrivial zeros of the 
Riemann zeta function have the form ^ + ia n , lying on a line pQ. This 
beautiful statement got related to mechanics by the conjecture of Hilbert 
and Polya, as a result, a search is on for a self- adjoint operator admitting 
real eigenvalues {cr n }. Perhaps the most striking work in this direction is due 
to Connes [2] who constructed a classical eigenvalue problem with a Perron- 
Frobenius operator and presented a spectral interpretation of Riemann zeros. 
In the realm of quantum mechanics, observations on trace formula have 
given some important insights and several Hamiltonians have been discussed 

El El 

The connection of the RH with random matrix theory (RMT) is very 
deep. For the statistical description of level-sequences (or number-sequences) 
of nuclei, Wigner [7j introduced the subject wherein Hamiltonian matrices 
were constructed keeping in mind the underlying symmetries possessed by a 
physical system. Thus, an even-spin, time-reversal invariant system belongs 
to a Gaussian Orthogonal Ensemble (GOE) whereas a system violating time- 
reversal invariance (TRI) belongs to a Gaussian Unitary Ensemble (GUE). 
That the sequence {c n } actually has a spectral interpretation first came 
out of the seminal work by Montgomery [B] who found the two-point cor- 
relation of Riemann zeros and obtained exactly the result known for GUE. 
Since then, higher-order correlations among Riemann zeros have also been 
shown to correspond to GUE; within GUE, it is known that two-point cor- 
relation guarantees all the higher-order correlations as they are factorisable. 
Perhaps the most important single effort after the Montgomery's work is 
the marathon numerics by Odlyzko [2] who has decidedly shown that the 
distribution has the form exactly as in GUE. Due to these works, the Hamil- 
tonians being searched for RH are the ones where time-reversal invariance 
is broken [3J |5J E] . 

Let us focus on two main points on which all the works rest. Firstly, the 
reality of eigenvalues of a Hermitian operator and completeness of solutions 
of the ensuing eigenvalue problem would guarantee that the RH holds true. 
Secondly, due to the mathematical and numerical works on correlations, it 
is expected that the Hamiltonian underlying the RH breaks TRI. In this 
paper, we construct a pseudo-unitary ensemble of random matrices which 
has the spacing distribution exactly as in GUE. Since these systems usually 
correspond to the physical situation where TRI and parity are not individ- 
ually preserved, the finding presented below suggests that the Hamiltonian 
underlying RH could also be pseudo-Hermitian. After this demonstration, 
we shall provide further reasons that attest to the above statement. 

A Hamiltonian H is called pseudo-Hermitian [141 115j if rfHrj 1 = for 
some metric rj. If E m and E n are two eigenvalues of H, it is known that JS] 

(E* m -E n ) (V* m \r,* n )=0, (1) 
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implying that if eigenvalues are real and different, the eigenstates are or- 
thogonal as (\P£j77\I/ n ) =G n <5 m>n . If an eigenvalue is complex it will have 
a zero pseudo-norm as = (^* |t/^ n ) = 0. The vanishing of pseudo-norm 
means that the eigenvector is null. If a Hamiltonian, H is symmetric under 
a joint action of parity P : x — > —x and time-reversal T : i — > —i, i.e, (PT) 
H (PT)" 1 =H [IH] then we have real eigenvalues if the eigenstates, ^ n are 
also the eigenstates of PT, otherwise the eigenvalues are complex conjugate 
pairs. For PT-symmetric Hamiltonians we have 

(E* m -E n ) (** T \* n )=0. (2) 

When eigenvalues are real and distinct, the eigenstates are orthogonal as 
(*£ T |^ n ) =£ n 5 m , n . 

Remarkably, PT-symmetric Hamiltonians are found to be pseudo-Hermitian 
: P H P _1 = Ht. Pseudo-Hermiticity has been recast |15| in terms of PT- 
symmetry. Given a pseudo-Hermitian Hamiltonian, one can construct gen- 
eralized P and T. 

The operator D= e*^ is pseudo-unitary in accordance with 

D f =f/D" 1 Tr 1 . (3) 

The eigenvalues of D are either on the unit circle or of the type : | A1A2I = 1- 
It is only recently that a random matrix theory has been presented for a 
statistical study of pseudo-Hermitian Hamiltonians ^2] . Only 2x2 matrices 
have been studied. To summarize briefly, two cases are found : one with a 
linear (with more slope than that of goe) level repulsion and the other where 
as the spacing s becomes small, level-spacing distribution ~ slog i |T2]. We 
call these ensembles as Gaussian pseudo-orthogonal ensemble (GPOE) 
and Gaussian pseudo- unitary ensemble (GPUE), respectively. The essence 
of these two results is that they show much weaker level-repulsion at small 
spacings than those of the known ensembles of Wigner and Dyson. 
Let us now consider the Hamiltonian matrix 



H = {H«} 



a + b (c + id)/e 
(c — id)e a — b 



,b,c,d being real. This is pseudo-Hermitian with respect to a metric 



n 



e 
l/e 



(4) 



(5) 



which gives rise to a positive definite pseudo-norm (1) It is due to this 
property such Hamiltonians as (4) are called quasi-Hermitian (see Scholtz 
et al. in |14j). The eigenvalues of H are given by 



E± = a ± Vb 2 + c 2 + d 2 . (6) 
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Consider that the matrix H is drawn from an ensemble of random matrices 
with a Gaussian distribution given by [7j 



P(H) 



, _ i tr H f H 

Ne 2^ 



Accordingly, the joint probability distribution of a,b,c,d is 

(c 2 + d 2 " 



P(a, b, c, d) ~ exp 



We know that three-parameter unitary matrix, U(0, 0, ■0) 



U 



cos 6» -sinfle** 
sin 0e~*^ e _M ^ cos # 



(7) 



(8) 



(9) 



constitutes a Lie group. More importantly, the unitary matrix U can gener- 
ate all the Hermitian 2x2 matrices of the general type (e = 1, in (5)) with 
any arbitrary value (including zero) of ip. Only two continuous parameters 
(9, 4>) suffice for this purpose. Inspired by this, we construct the following 
matrix D: 



cos V 
sin0e _i *e 



cos 9 



(10) 



which is pseudo-unitary with respect to rj. As per our design D matrix would 
generate all possible H of the type (4) as 



D diag(E + ,£;_)D~ 1 = H, 

which gives us the following relations : 

_£/_)_ + E_ _£/_)_ — E_ 



a 



2 ' 2 
c = sin 29 cos (j), d 



cos 29, 
E+ - E- 



sin 29 sin < 



, • - 2 (12) 

Writing e = e~ 7 , and calling t = E + + E_, s = E + — EL, we have P(a, b, c) 
going over to 



P 7 (s, t, 9, 4>) ~ exp 



- —^7 — —7 cos 2$ 7T cosh 27 cos sin 1 

4cr 2 4cr 2 a 2 



(13) 



via a Jacobian, J7 = g2 s ^ n 20 . Next, integrating over i, 6* and (f>, we have the 
un-normalised nearest-neighbour spacing distribution given by 



-P 7 (s) ~ s exp 



p 2 s 2 



Erfi 



V2cj 



(14) 



where 



Erfi(x) 



y dye x2y \p = v / cosh2 7 /2,g = v / cosh2 7 - 1/2 (15) 
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Figure 1: Plot of P-y{x) (solid line) (16) for three values of 7. Dashed line 
denotes Pgue(x). 



We write now the normalized nearest-neighbour spacing distribution in terms 
of a dimensionless variable, x = -Ar where (s) is the mean level spacing. With 



2 r tanh-^q/pU 
" 7 V 4pq J 



= ^cosh2 7 x e _ p2a>2 



Ern(a 7 gx) 



(16) 



Note that the limiting value of the square-bracketted term is and ccq = 
— y=. For an arbitrarily small 7, the matrix H is pseudo-Hermitian. Actually, 
even for 7 as much as 1/2, the difference between -P 7 (x) and Pqu'e(x) = 

||x 2 e — is hardly appreciable (see Fig. 1(a)). 

Returning to the discussion of the RH, with this example-ensemble, the 
scope of search of the Hamiltonian for the RH widens. The Hamiltonian 
in question relevant for the RH could be pseudo-Hermitian. Although the 
example presented is non-generic, so could the Riemann Hamiltonian be. 
There is nothing that suggests generic nature of the Hamiltonian, particu- 
larly in the light of our illustrative example. 

Our suggestion is also well supported by the spectral interpretation of 
the Riemann zeros ensuing from Connes' work According to Connes, 
the zeros form an absorption spectrum in the sense that the wavefunctions 
corresponding to the eigenvalues a n is "zero" . We know that the eigenvalues 
of a pseudo-Hermitian operator are either real or complex-conjugate pairs. 
Thus, we suggest the possibility of the Riemann zeros ± ia n } to be the 
complex- conjugate-pair eigenvalues of an unknown pseudo-Hermitian oper- 
ator where it would be automatically guaranteed that the eigenvectors are 
null. Our paper suggests this central message. 
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Let us demonstrate our point heuristically by taking a simplistic and 
trivially PT-symmetric Hamiltonian as 

Hpj = —ix p. (17) 

The eigenvalues like \ ±it n will be supported with ^ n (x) = Nx~^ ltn such 
that ^ n (±oo) = and Hpj^„, = ± it n )^ n (x). Very importantly notice 
that the eigenvalues are complex conjugate pairs and the eigenfunctions of 
Hp-p are not the simultaneous eigenstates of the antilinear operator, PT. As 
stated earlier, this situation is referred to as spontaneous breaking of PT- 
symmetry [TO]. Check that PT$> n (x) = ^* n (-x) ^ c^ n (x). Next, whether 
^ + it n are bona fide discrete eigenvalues and whether t n would coincide 
with a n (RZs) are of course the most crucial questions. 

Our simple Hamiltonian in (17) mimics the Hamiltonian of Berry and 
Keating [3| 

H B K = *P-^ (18) 

which, in turn, has been inspired by the work of Connes. This is Hermitian 
and it also breaks time-reversal symmetry. Berry and Keating [3J |lj have 
studied the the semiclassical trace formula for their Hamiltonain (18) vis-a- 
vis very interesting properties of C(z) and reported a shortcoming of (18) in 
this regard. Also they speculated (21 that the Hamiltonian (18) along with 
and extraordinary boundary condition on the wavefunction would yield ±<7 n 
as eigenvalues. This boundary condition is unfortunately not known so far. 
Also, if a n is an eigenvalue, apparantly there is nothing to ensure that — a n 
would also be an eigenvalue. 

When we diagonalize the matrix for Hg^ using the one-dimensional 
Harmonic Oscillator basis by using the creation and innhallation operators 
as : x=(a+a^)/\/2 and p=i(a^-a)/ \/2, we find that eigenvalues very crucially 
depend upon the size of the basis (say iV)! This is how we conclude that 
HbK does not even possess a discrete spectrum. So is the fate of our toy 
model Hp-p (17), this however is only heuristic. These simple findings are 
for the most ordinary boundary condition where the eigenfunctions vanish 
at ±oo. The real part turns out to be 1/2, and this would change as soon 
as the boundary conditions are disturbed. 

The classical analogue of the Hamiltonian (18) is known to be scaling 
type (as x— > K x, p— > p/K), therefore the complex scaling of co-ordinate 
can not be employed to study its resonances. However, its canonically- 
transformed Hamiltonian, H = (p 2 — x 2 )/2 is very well- studied for its 
resonances and these are well-known as ±i(n + ^) (with h = 1) not showing 
any connection with a n . 

Okubo has considered the Hamiltonian, 

i 

H Okubo = -PxPy - (1 - /3)xPx - fiyPy + ^ (19) 
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that is both Hermitian and time-reversal breaking. It is in two-dimensional 
Euclidean space with boundary conditions on the eigenfunctions : 

H Okubo^( x > V) = A ^(^ y); i/>(x, 0) = 0, (20) 

and ip(x,y) rapidly decreasing at infinity. We have again constructed the 
Hamiltonian matrix for (19) in the harmonic oscillator basis and diagonalised 
the matrices to find the eigenvalues. In x it is a H.O. and in y it is half- 
H.O. model as per We find that the eigenvalues are not stable with 
the size of the matrices, thus indicating that there is no discrete spectrum 
supported by the Hamiltonian. However, the possibility of the Riemann 
zeros to correspond to resonances remains with this Hamiltonian. It may 
again be non-trivial to investigate its resonances. 

The Hamiltonian suggested by Castro et al. ^7] is like H^qm = ix p + 
g(x), where g(x) is a special function and x > 0. Once again, by finding the 
matrix elements employing half H.O. basis, we do not find discrete specrum 
as the eigenvalues keep changing with the size, N, of the basis. 

We have found that the popular Hamiltonians in the context of the RH do 
not even possess a discrete spectrum. From the random matrix ensemble of 
pseudo-Hermitian matrices presented here exhibiting GUE statistics and by 
a re-interpretation of Connes' work, we have suggested that the Hamiltonian 
relevant to the RH could be pseudo-Hermitian. 
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